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Abstract 

We develop the theory of fermion induced phononlike excitations of the instanton liquid. It 
suggests a mechanism of current quark mass generation which is easily understandable by cal- 
i—j | culating the corresponding functional integral in the tadpole approximation. We systematically 

\^ • study the quark condensate excitations influenced by the phononlike modes and rederive the 

CN ■ relation Gell-Mann-Oakes-Renner with realistic pion mass. The picture of u-meson as being 

mixed with the soft scalar glueball-like excitation is discussed. 
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*q} I. Introduction 
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Nowadays, there are no doubts the model of QCD vacuum as the instanton liquid (IL) is the most 
practical instrument on the chiral scale of QCD. It provides, as the lattice calculations recently con- 
firmed, not only the theoretical background for describing spontaneous chiral symmetry breaking 
(SCSB) but is mostly powerful in the phenomenology of the QCD vacuum and in the physics of light 
quarks while considered to propagate by zero modes arising from instantons. The origin of gluon and 
chiral condensates turns out in this picture easily understandable and both are quantitatively calcu- 
lated getting very realistic values defined by Aqcd and parameters of instanton and anti-instanton 
ensemble, for example, — % (ip'ip) ~ —(250 MeV) 3 . Moreover, the scale for dynamical quark masses, 
M ~ 350 MeV, naturally appears and pion decay constant, f n ~ 100 MeV, is then transparently 
calculated. 

Another significant advantage of this approach is that the initial formulation starts basically from 
the first principles and subsequent approximations being well grounded and reliably controlled are 
plugged in [|IJ, 0. It becomes clear especially in recent years when the impressive progress has been 
reached in understanding the instanton physics on the lattice |§ . Further we are summarizing several 
things we have learned thinking of the IL theory and trying to answer the challenging questions. 

Let us start on that stage of the IL approach when its generating functional has already been 
taken as factorized one into two factors 

Z = Zg ■ Z^j , 

where eventually Z g provides nontrivial gluon condensate and the fermion part Z q is responsible to 
describe the chiral condensate in instanton medium and its excitations. It is usually supposed the 
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functional integral of Z g is saturated by the superposition of the pseudo-particle (PP) fields which 
are the Euclidean solutions of the Yang-Mills equations called the (anti-)instantons 



A' 
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Here A^x; 7A denotes the field of a single (anti-)instanton in singular gauge with 4iV c (for the SU (N c ) 
group) coordinates 7 = (p, z, U) of size p centred at the coordinate z and colour orientation defined 
by the matrix U. The nontrivial bloc of corresponding N c x N c matrices of PP is a part of potential 

^(2; 7) = f ^ L %^- i — 2 lPr a U, y = x - z , a = 1,2,3 , (2) 

9 y y +p 

where r a are the Pauli matrices, rj is the 't Hooft symbol ||], g is the coupling constant and for anti- 
instanton 77 — 77. For the sake of simplicity we do not introduce the distinct symbols for instanton 
(N + ) and anti-instanton (iV_) and consider topologically neutral IL with N + = iV_ = N/2. Utilizing 
the variational principle the following estimate of Z g was found @ 

Z g ~ e~ {s) 

with the action of IL defined by the following additive functional [] 

(S) = dz dp n(p) s(p) . (3) 



The integration should be performed over the IL volume V along with averaging the action per one 
instanton 

s(p) = (3(p) + 5 ln(Ap) - lnp 2N ° + ^p 2 J dp, n{ Pl )pl (4) 
weighted with instanton size distribution function 

n(p) = C e~ s{p) = C p - 5 ^ e -/3(p)-P 2 /7 ; (5) 

6-4 , llA^ c -2A^ f v 
z/ = ^— > & = o ' \P 
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where p 2 = J dp p 2 n(p)/n, n = J dp n(p) = N/V and Nf is the number of flavours. The constant 
C is defined by the variational maximum principle in the selfconsistent way and (3{p) = ^5- = 
— lnCAr c — 61n(Ap) (A = Ajjg = 0.92Apy) with constant CV C depending on the renormalization 

scheme, in particular, here C^ c ~ eX ^) j (N^2) ! ' ^ e P arame ^ ers P = P(p) an d P = (3+\iiCn c 

27 N 

are fixed at the characteristic scale p (an average instanton size). The constant £ 2 = -^2 c ^ 2 

characterizes, in a sense, the PP interaction and Eqs. ©,(§) and (|5]) describe the equilibrium state 
of IL. The minor modification of variational maximum principle (see Appendix) leads to the explicit 

form of the mean instanton size pA = exp j — ^7~rjj an d> therefore, to the direct definition of 

the IL parameters unlike the conventional variational principle which allows one to extract those 
parameters solving numerically the transcendental equation only. 

1 In fact, the additive property results from the supposed homogeneity of vacuum wave function in metric space. 
Eq. (||) looks like a formula of classical physics although it describes the ground state of quantum instanton ensemble. 
Intuitively clear, this definition will be still valid even when the wave function is nonhomogeneous with the nonuni- 
formity scale essentially exceeding average instanton size or, precisely speaking, being larger (or of the order) than 
average size of characteristic saturating field configuration. Then each instanton liquid element of such a distinctive 
size will provide a partial contribution depending on the current state of IL, see next Section. 



The quark fields are considered to be influenced by the certain stochastic ensemble of PPs, Eq. 
flU), while calculating the quark determinant 

Z f ~ / D^Lhj) « eSWA )) A . 



Besides, dealing with dilute IL (small characteristic packing fraction parameter n p 4 ) one neglects 
the correlations between PPs and utilizes the approximation of N c — > oo where the planar diagrams 
only survive. In addition the fermion field action is approached by the zero modes what means the 
quark Green function is considered as the superposition of the free Green function So = (— id)^ 1 and 
fermion zero modes — z) which are the solutions of the Dirac equation i(D(A±) + m) &± = 
in the field of (anti-)instanton A± centred at z, i.e. 

$ ± $^ 

S± = S ; , (6) 

i m 

here m is the current quark mass and the quark zero mode possesses the following analytic form 

l±7s- 
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with the colour c, d and the Lorentz i,j indices and the antisymmetric tensor e. In fact, there exists 
the singular term in the limit m — ► but it is selconsistently fixed by the saddle point calculation of 
quark determinant Z^. In spite of the fact the exact Green function (at m^O) including the terms 
of the whole series is well known [0 and, moreover, the Green function of instanton molecule has 
been also established |J, the simple zero mode approximation (|) is still the most practical in the 
concrete evaluations. In particular, at Nf = 1 the quark determinant reads 

Z^ ~ fDtfDif, e w ydx^(x)id*l>(x)} (J^j + , (7) 

Y ± — i dz dU dp n(p)/n / dxdy ip\x) id x $±{x — z) $±(y — z) id y ip{y) , 



where the factor M makes the result dimensionless and is also fixed by the saddle point calculation. 
Amazingly this relatively crude approximation turns out so fruitful to develop (even quantitatively!) 
the low energy phenomenology of light quarks. The generating functional beyond the chiral limit 
was obtained in Ref. 0. 

Thus, the IL approach at the present stage of its development looks very indicative, well the- 
oretically grounded and reasonably adjusted phenomeno logically. The proper form of generating 
functional obtained and its verisimilar parameter dependence indicated provide enough predictive 
power and justify, hence, the approximations made. It dictates the improvements to be done within 
the approach but, on the other hand, calculating some corrections has clearly no serious prospect. 
Taking this message as a guiding one we are going to demonstrate that amplifying the approach with 
an inverse influence of quarks upon the instanton ensemble which is intuitively small effect leads, 
however, to rather unobvious important conclusions. 

We describe this influence (not getting beyond IL and SCSB approach) as a small variation of 
instanton liquid parameters 5n and 8p around their equilibrium values of n and p being in full analogy 
with the description of chiral condensate excitations. Indeed, the result of nontrivial calculation of 
the functional integral (treating substantially from physical view point the zero quark modes in the 
fermion determinant) comes to 'encoding' the IL state with just those two parameters. Moreover, 
the IL density appears in the approach via the packing fraction parameter np 4 only (clear from 
dimensional analysis) what means one independent parameter existing in practice. It is just the 
instanton size. The analysis of the quark and IL interaction is addressed in this paper developing 
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our idea || of phononlike excitations of IL resulting from the adiabatic changes of the instanton 
size. Thus, we will suppose the inverse quark influence should be described by these deformable 
(anti-) instanton configurations which are the field configurations Eq. @ characterized by the size p 
depending on x and z, i.e. p — > p(x, z). 

The paper is organized as follows: in Section II we discuss the modification of quark determinant 
when the functional integral is saturated by the deformable modes at the minimal number of flavours. 
Then in Section III we develop the approximate calculation (tadpole approximation) which is based 
on the saddle point method and the corresponding iteration procedure. Section IV is devoted to 
the generalization for the multiflavour picture. The meson excitations of quark condensate and 
calculation of the Gell-Mann-Oakes-Renner relation when it is provided by the mechanism of quark 
current mass generation related to the phononlike excitations of IL are analyzed in Section V. The 
paper includes also Appendix where the fault finding reader gets a chance to control the explicit 
formulae of the IL parameters and to improve our calculations if is able. 



II. Supplementing phononlike excitations 

Apparently, the gist of what we discuss here could be illuminated in the following way. Saddle point 
method calculation of the functional integral implies the treatment of the action extremals which are 
the solutions of classical field equations. For the case in hands the action S[A, ip^,i/>] is constructed 
including the gluon fields A^, (anti-)quarks ip',ip and extremals, which are given by the solutions 
of the consistent system of the Yang-Mills and Dirac equations. As a trial configuration in the IL 
theory the superposition of (anti-)instantons which is the approximate solution of the Yang-Mills 
equations (with no reverse influence of the quark fields) and an external field for the Dirac equation 
simultaneously is considered. We believe it is reasonable to utilize the deformable (crampled) (anti- 
)instantons A±(x\ j(x)) as the saturating configurations. They just admit of varying the parameters 
j(x) of the Yang-Mills sector of the initial consistent system in order to describe the influence of 
quark fields in the appropriate variables for the quark determinant. 

Taking the action in the form S[A±(x, , y(x)),iJ)',il)] we would receive the corresponding varia- 
tional equation for the deformation field j(x) which approaches most optimally (as to the action 
extremum) PP at nonzero quark fields. In the field theory, for example, the monopole scattering 



@, the Abrikosov vortices scattering |10[ are treated in a similar way. However, for IL we avoid 
the difficulties which come with solving the variational equations if we consider the long-length wave 
excitations only with the wave length A much larger the characteristic instanton size p. Indeed, it 
looks possible because we are searching the kinetic energy of the deformation fields (the one particle 
contributions) and consider the pair interaction which develops the contact interaction form being 
calculated in the adiabatic regime 

Let us remind first of all that deriving Eq. ([3]) we should average over the instanton positions in 
a metric space. Clearly, the characteristic size of the domain L which has to be taken into account 
should exceed the mean instanton size p. But at the same time it should not be too large because 
the far ranged elements of IL are not 'causally' dependent. The ensemble wave function is expected 
to be homogeneous (every PP contributes to the functional integral being weighted with a factor 
proportional to ~ 1/V, V = L A ) on this scale. The characteristic configuration which saturates 
the functional integral is taken as the superposition Eq. ([]]) with N of PP in the volume V. It is 
easy to understand that because of an additivity of the functional Eq. (|]) describes properly even 
non-equilibrium states of IL when the distribution function n(p) does not coincide with the vacuum 
one and, moreover, allows us to generalize it for the non-homogeneous liquid when the size of the 
homogeneity obeys the obvious requirement A > L > p. Besides, we should deal with the weak 



2 Then we are allowed to take the slowly changing deformation field beyond the integral while calculating the action 
of deformed instanton. 
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(comparing to the instanton G^ u ) fields g^ u (g^ u <C G^ u ) and long- length wave (on the scale of 

dp(x, z) 



the instanton size p) perturbations Q% — ^ 0(1) when the ensemble saturating the functional 

integral is very close to the vacuum (instanton) one. Then this smoothness or the adiabatic change 
of instanton size practically dictates another essential simplification defining everywhere the field in 

OP\X Z) op(x z 
the center of instanton V, ' — - ~ \ ' — 



characteristic deformation field and the exact 

x=z 



correction to the potential 



y 



dp 
dx v 



where $^(x, z) = ^£S2^1 — ^_ Keeping ourselves within the precision accepted here we 

could take p(x,z) ~ p, i.e. $^(x, z) = &u V (x — z). On the other hand the most general form of 
the correction to the instanton field influenced by the quarks might be calculated if the gluon Green 
function in the instanton medium is known 

a*(x,z) = Jd£ Df u (x -z^-z) J b u ^ - H ) , (9) 

where J h v is the current of external quark source, z$ belongs to the region of long-length wave 
disturbation and, at last, D°f v is the Green function of PP in the instanton medium. In fact, this 
function is not well defined but, seems, for the case in hands we could develop the selfconsistent 
way to calculate the regularized Green function. The nonsingular propagator behaviour in the soft 
momentum region is defined by the mass gap of phononlike excitations. Fortunately, the exact 
form of the Green function occurs unessential here (we are planning to return to the problem of 
regularized Green function calculation in the forthcoming publication). In the coordinate space it is 
peaked around the averaged PP size being in nonperturbative regime and, hence, integral Eq. @ is 
appraised to be 

a*(x,z)~F£(x-z) Jt{z- H ) . (10) 
Appealing now to Eq. <ffl) we are capable to get immediately for p v the following equation 



%te - *) = f;1{x - z) 4( z - ^) . 



In fact, the current J might be taken constant in the long- length wave approximation. Then 



dp(x, z) 



dx 



dp(z) 



dz 




neglecting the gradients the following change appears to be justified everywhere 

since there is no other fields in the problem at all (in the adiabatic approximation). The deformable 
mode contribution to the functional integral (when the corrections coming from the deformation 
fields of PP are absorbed) may be estimated as 

(S) - J dz J dp n{p) 

where k is the kinetic coefficient being derived within the quasiclassical approach. Our estimate of 
it gives the value of a few instanton actions k ~ c j3 with the coefficient c ~ 1.5 4- 6 depending 
quantitatively on the ansatz supposed for the saturating configurations. Although this estimate 
is not much meaningful because there is no the vital k dependence eventually (becomes shortly 
clear). Thus, this coefficient should be fixed on a characteristic scale, for example k ~ n(p) if we 
are planning not to be beyond the precision peculiar to the approach. Actually, it means adding 
the small contribution of kinetic energy type to the action per one instanton only. Such a term 



results from the scalar field of deformations and affects negligibly the pre-exponential factors of the 
functional integral. In one's turn pre-exponential factors do the negligible influence on the kinetic 
term as well |§. If we strive for to be within the approximation we should retain the small terms 

ds(p) „ 
— = supposing 



of the second order in deviation from the point of action minimum only 
approximately 

s(p) ~ s(p) + , 

= = and the scalar field <p = 5p 
p 2 

deviations from the equilibrium value of p c = p ( 1 — ^7 ) — 
field is described by the following Lagrangian density 



P=Pc 



where s <y2 \p) 



d 2 s{p) 



P ~ Pc 



p — p is the field of 



p. Consequently, the deformation 




£ = 

with the mass gap of the phononlike excitations 

sW(p) 



2,„2 



M 



K 



+ M z <p 



up 2 



which is estimated for IL with N r 



3, for example, in the quenched approximation to be 
M fa 1.21 A 



if c = 4, p A fa 0.37, (3 fa 17.5, n A" 4 fa 0.44 (for the details see the tables of Appendix). The 
deformation fields (p (with corresponding Jacobian) contribute to the generating functional on the 
same footing as the quark fields and it looks like 



Z ~ 

9 



Dip 



6A 



5ip 



exp 



dz 




+ mV 



in full analogy with the fields ip\ ip entering with the functional measure Dip\ Dip. Actually the 
Jacobian contribution [] should be omitted in what follows as discussed above. 

Analyzing the modifications arising now in the quark determinant we take into account the 
variation of fermion zero modes resulting from the instanton size perturbed 



where ^± (u, p c 



®±(x - z,p) ~ $±(x - z,p c ) + $i 1} (x - z,p c ) 5p(x,z) , 

and because of the adiabaticity it is valid Sp(x, z) 



d$±(u } p) 
cfp 



8p(z,z) 



p=Pc 



(p(z). The additional contributions of scalar fields generate the corresponding corrections in the 
factors of the kernels Y ± of Eq. (|7|) which are treated in the linear approximation in (p, i.e. 



id x $±(x - z, p) <& ] ±(y - z, p) id y 
here we introduced the notations 



r±(x,y,z,p c ) + T ( l\x,y,z,p c ) 



<p(z) 



(12) 



±(x,y,z,p c ) 



id x $±{x - z, p c )^±{y - z, p c )id. 



V J 



3 Generally, the deformation field p v and integration variable (|S|) are related via the rotation matrix: il a b$* !/ (a; — 
z)p v = aJJ and in the long- length wave approximation $ might be constant <I>^„(0) (x ~ z). With the rotation matrix 
spanning the colour field — ri _1 a M on the fixed axis (on the z axis for SU(2) group, for instance) we can conclude 
that the vectors and p v are, in fact, in one to one correspondence (of course, being within one loop approximation 
and up to this unessential colour rotation). Thus, the Jacobian occurs an unessential constant. 



r^\x,y,z,p c ) = id x $ ( ±'(x - z,p c )&±(y - z,p c ) id y + id x $±(x - z,p c ) &±'(y - z,p c ) id y 

(the gradients of scalar field <p are negligible according to the adiabaticity assumption again). It is a 
simple matter to verify that the right hand side of Eq. (12), being integrated over dzdU, generates 
the following kernel (in the momentum space) 



t(i), 



Nr. 



(2vr) 4 6{k-l) l0 {k,k)+ 7l {k,l) p(k-l) 



(13) 



with j Q (k,k) = G 2 {k) , G(k) = 2n Pc F{kp c /2) , ^(k,l) = G{k)G'{l) + G'{k)G{l) , 

, F(x) = 2x [I (x)K 1 (x) - h(x)K (x)} - 2 I 1 (x)K 1 (x), where K { (i = 0, 1) 



G'(k) 



dp 



P=Pc 



are the modified Bessel functions. 

In fact, the functional integral of Eq. (|7|) including the phononlike component may be exponenti- 
ated in the momentum space [] with the auxiliary integration over the A-parameter (see, for example 

E) 



Z,^J- expjiVln 
x / Dip^Dip exp 



N 



iXVM 
dkdl 



N\ x 



%p\k) 



iX 



(2nf5(k -l)[-k + — 7o (fc, k) + — 7l (fc, I) p(k - I) 



N 



Nr. 



(we dropped out the factor normalizing to the free Lagrangian everywhere). It is pertinent to mention 
here the Diakonov-Petrov result comes to the play precisely if the scalar field is switched off. 

In order to avoid a lot of the needless coefficients in the further formulae we introduce the 
dimensionless variables (momenta, masses and vertices) 



— — > k 
2 



Mp c 



M , 70 -> Pc7o , 7i 



Pell 



the fields in turn 



v{k ) _> (nK )-l/2 p ^ {k) j m _> p l/2 m 



(14) 



(15) 



and eventually for A we are using p = Qjy- Then the generating functional takes the following form 

dk 1 



Z ~ J dfiZ'g J D^Dip Dip exp |-iVln/x - J 



7T 4 2 



4 + M ) <p(k) x 



x exp 



dkdl 



71 



^(k) 2 



7T 4 5(k - l)(—k + ipjoik, k)) + 



ip 



mli(kj) p(k-l) 



(16) 



where Z" is a part of gluon component of the generating functional which survives after expanding 
the action per one instanton Eq. (|IID . The functional obtained describes the IL state influenced 
by the quarks when all the terms containing the scalar field are collected (see also Appendix). As 

In the metric space we have the nonlocal Lagrangian of the phononlike deformations f{z) interacting with the 
quark fields ijv , ip, i.e. 



C = 



dx r/y (x)id x ip(x) — dz 



dip 

dz 



+ M 2 <p 2 (z) } + 



dxdydz dU ip f (x){T±(x, y, z, p c ) + (x, y, z, p c ) (p(z)}ip(y) . 



The physical meaning of the basic phenomenon behind this Lagrangian seems pretty transparent. The propagation 
of quark fields through the instanton medium is accompanied by the IL disturbance (the analogy with well known 
polaron problem embarrasses us strongly in this point). 
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mentioned above we believe this influence analogous to the reversal impact of phononlike deformations 
on the quark determinant does not considerably change the numerical results of the IL and SCSB 
theory. But it is invisible from Eq. (jlBf) directly how this smallness is reasoned. The scalar field 
enters the generating functional formally at the same order of the /i-expansion as the term providing 
SCSB (compare the second and third terms of the second exponential function in Eq. fll6|) ). The 
suppression arises because the dominant contribution of the Yukawa interaction comes from the 
quark field condensate which maintains the additional fi smallness. This result prompts, in fact, 
rather natural scheme of the approximate calculation of the generating functional 

If we believe, for example, the first term ascribes some gluon component of the generating functional 
then being linearly dependent on the scalar field it produces the small shift from the equilibrium 
value of the instanton size (p c ~ p) (p — tp' + Sip. And the shift 5<p generates the mass term in the 
quark sector what means the scheme should be pushed forward to the following expression [] 

^ty y _ (y _|_ Sip) + _ (^t^)} {^pf + Sip) . 

Thus, we face the conventional mechanism of the mass generation when it is related to the insignificant 
variation of the equilibrium instanton size p c (or, in usual terms, to falling the scalar field condensate 
down) produced by the quark condensate. Moreover, it is clear if the variations happen to be of 
the opposite sign then the field 75^ (the chiral partner of the field ip) should develop the true 
mass. Before turning to the mechanism of the mass generation below let us calculate the quark 
determinant integrating formally over the scalar field Dip, then find the quark Green function in the 
tadpole approximation and formulate the equation for the saddle point. It is supposed the phonon 
component contribution does not affect substantially the results of the SCBS theory. Here we are 
interested in projecting upon the scale inherent for the scalar field. 

III. Tadpole approximation 

The integration leads us to the four fermion interaction and the functional integral can not be 
calculated exactly. However, due to smallness of scalar field corrections we may find the effective 
Lagrangian substituting the condensate value in lieu of one of the pairs of quark lines (see Fig. 1.) 




fermion fields s/>, ^ 
phononlike field ip 



Figure 1: The tadpole contribution. 

^{k)ip{l) -> (ip\k)ip{l)) = -7r 4 5{k - I) Tr S{k) . 

In such an approach the diagram with four fermion lines in the lowest order of the perturbation 
theory in /z is reduced to the two-legs diagram with one tadpole contribution (there are two such 

5 Let us emphasize the first term is responsible, in a sense, for the quark condensate fluctuations what allows the it 
meson to gain its mass. In the second term such fluctuations are suppressed. 



Q 



contributions because of two possible ways of pairing) 



-4 

np 



\ (ip)_ r dkdl dk'dl' 7i(M)7i(fc , r) ^ {k)m ^ {k > )m (^(A; - i) ip(k' - I')) ~ 

K J 7T 



' ^ / ^ 7l (fc,fc) ^(fc)^(fc) / 4 7i(M) Tr D(0) , 

J 71 J 7T 



np 4 n 



where the natural pairing definition was introduced 

{<p(k)<p(l))=**6(k-l)D{k), D{k) = 4 I . 
It is obvious the factors surrounding ijfi \k)ip{k) has a meaning of quark mass 

m ' w = -"<*•*> w$n!'^^ Tr s « • (17 > 

(the initial mass term contains the factor 2 when the dimensionless variables are utilized, i.e. takes 
a form 2im/.) We are treating it as the current quark mass since the further calculations show its 
magnitude gets just within the scale interval commonly accepted for the current quark mass and, 
moreover, it appears in the Gell-Mann-Oakes-Renner relation. 

The contribution of the graph with all the quark lines paired (see Fig. 2) should be taken into 





Figure 2: Second order contribution, 
account at the same order of the p expansion while calculating the saddle point equation 
2p 2 



np A n 



dk 

— 7i(M) Tr S(k) 

7T 



1 2 



7r 4 <S(0)D(0) 



1 p 2 V_ 

o -4 -4 

2 np v p 



dk 

- 7 i(M) Tr S(k) 

7T 



1 2 



-KYx in the dimensionless units. 

71 p 

Then the quark determinant after integrating over the scalar fields reads 



Here we used the natural regularization of the 5-function 5(0) 



dp exp < — N In p + . 

{ np v p 



/r j o JV 2 V t dk 

dp / D^Di/j exp \-N\np + -j p A c 2 {p) + / ^(k) 2 [-k + iF(k)] ip{k) 

J [ np v p J ix 

2N 2 V 4 2 , 



V r dk 

- //' <*(,,)+ \ —,Tt \n[-k + ir(k)} 



where the vertex function is defined as 

F(k) = p 7o(fc, k) + m f (k) , 
and we introduced the function c(p) convenient for the practical calculations 

<*0 = -^ /*•«(*,*> i>5<*). 



Figure 3: The function c(/i) at Nf — 1. 



As it is clear from Eq. (|T^) the Green function of the quark field is self-consistently defined by the 
following equation 

2 [-k + iT(k)) S(k) — -1 . 

Searching the solution in the form 



S{k) = A(k) k + i B{k) 

we get 



i r(*) 



2 P + r 2 (A;) ' w 2 A; 2 + T 2 (A;) ' 

Using Eq. ([H]) and the definitions of T(k) and B(k) we have the complete integral equation 

T(k)=^ l0 (k,k) + -^-^ 2 7l (k,k) 1(0 



npV " ,iV ' ' J 7T 4 ,±v ' v z 2 + r 2 (/) 

which drives to have the convenient representation of the solution 

Nr. 



T(k) = fj, 7o(fc, fc) + — ^ /i 3 c(/i) 7i(/c, fc) . 



What concerns the function c(/i) it is not a great deal to obtain 

n i r dk r(k) 



/i 7 tt 4 m ' v A; 2 + r 2 (A;) ' 

and, therefore, the complete integral equation for the function c(/i) which is shown in Fig. 3 for 
Nf = 1. Let us underline the iV/- dependence of the c(/i) function in the interval of /i determined by 



i n 



Figure 4: The branches of the function c(p). 



saddle point value is unessential. Then we easily obtain for the current quark mass 

N u 3 

m f (k) = Z- c(/i) 7 i(M) , (19) 
np v 

and see the cancellation of the kinetic coefficient k in m/. Thus, it means the precise value of the 
coefficient is unessential as declared. 

We have the following equation for the saddle point of the functional of Eq. (|T8|) 



2Nc fdk M F 2 (A;)]; | 2N* p, r _, 2 
np 4 J 7r 4 A; 2 + r 2 (fc) np 4 z/ np 4 



+ -4^ ^4 \f <?Wu = 1 , ( 2 0) 



where the prime is attributed to the differentiation in p. 

It results from assuming the stationary IL parameters, what is not accurate. We should include 
another effect produced by the shift of equilibrium instanton size thanks to the quark condensate 
presence. The modification of the IL parameters {n(p), ■ ■ •) caused by the Yukawa interaction comes 
from simple tadpole graph of Fig. 5 in the leading order 



2 ip r dkdl , , i 
(np 4 Ky/ 2 J TT 8 



7 1 (ife, I) (-vr 4 ) 5(k - Tr S(k) <p(k - I) = A • y?(0) , 

(21) 



6 The solution is unique within the n interval of our interest. It is curious to notice when \x is larger than ~ 4- 10 2 
several solution branches for c(/i) emerge (see Fig. 4). Surely, it could be interesting to clarify if there is any 



correspondence between these branches and the saddle point equation resulting from Eq. (18). However, it is clear 
wittingly those objects should be heavier than the scale of several hundred MeV characteristic for SCSB. It is quite 
possible such solutions might be associated with some 'heavy' particles if they do exist. 



Figure 5: The tadpole. 



and this term generates the small shift of the equilibrium instanton size p c ~ p (remind here (p = p—p c , 
and (f(0) = J dz (p(z) is the scalar field in momentum representation) Q. The tadpole contribution 
could be absorbed in making more sophisticated variational procedure of saddle point calculation. 
Actually, it includes also the variation of the IL parameters as a function of p. But in practice it 
comes about highly effective to use the simple iterating procedure. At the first step this variation of 
the IL parameters is not taken into account and the saddle point p(p c ) is calculated from Eq. fl20|) . 
Then getting new IL parameters (see, Appendix) we have to resolve Eq. (|20"D again and etc. Five 
or six iterations are quite enough if we are satisfied with the same precision as in calculating the 
integrals. In the Table 1 the numerical results (M.S.Z.) are shown for Nf = 0, 1 (Nf = corresponds 
to the quenched approximation for the IL parameters) comparing to those of Diakonov and Petrov 
(D.P.) where the disturbance of instanton medium was not considered. 

Table 1. 



D.P. 


M.S.Z. 


N f p M(0) 


N f pt M(0) -ityty) m f 


5.68-KT 3 341 -(301) 3 

1 5.14-10" 3 376 -(356) 3 


5.51-10- 3 366 -(297) 3 4.29 

1 5.07-KT 3 385 -(328) 3 5.13 



The parameters indicated in the Table 1 are 
the dynamical quark mass 



M(0) = T(0) — 



[MeV] 



the quark condensate 



—i(ip>ip) = —i Tr S(x)\ 



2=0 



Nc 
4 



dk T(k) 
~^k 2 + T 2 (k) 



[MeVf 



and, at last, the current quark mass mf [MeV], which is defined now being mixed with the quark 
condensate as 

[dk m f (k) Tr S(k) 
mf ~ Jdk Tr S(k) ■ (22) 

It is just what is dictated by the Gell-Mann-Oakes-Renner relation which we calculate below. 
Through this paper the value of renormalization constant is fixed by A = 280 MeV. Then the IL 
parameters are slightly different from their conventional values p ~ (600 MeV)" 1 , R ~ (200 MeV)" 1 
(see the corresponding tables in Appendix). However, with the minor A variation the parameters 
could be optimally fitted. As expected the change of quark condensate is insignificant, the order of 



7 This shift in phononlike component of Lagrangian —4 M^- ip 2 + A 



_ 4 M 2 ,J2 



if — Sip and Sip = - generates the mass term in quark sector m^(k) = 



■ ■ with the definitions 
i \Tpi 7l(k>k) ^> wmcn has 



obviously the same form as that of Eq. 



(np k) 



several MeV, what hints the existence of new soft energy scale established by the disturbance which 
accompanies the quark propagation through the instanton medium. 



IV. Multiflavour approach 

In order to match the approach developed to phenomenological estimates we need the generalization 
for Nf > 1. Then the quark determinant becomes P] 



Z i> - J D^Dip exp I J dx X>}(x) idiP f (x)\ [yM^) 



N- 



N f 

Y± = i N f / dz dU dp n(p)/n J[ / dx f dy f ipj(x f ) id Xf $±(x f - z) $±(y/ - z) id Vf ipf(y f ) . 

/=i 

With phononlike component included every pair of the zero modes ~ $ $t acquires the additional 
term similar to Eq. (|T2|). The appropriate transformation driving the factors Y to their determinant 
forms H is still valid here since the correction term differs from the basic one with the scalar field 
p. The complete integration over dz leads (in the adiabatic approximation p>(x,z) — > p>(z)) to the 
transparent Lagrangian form with the momentum conservation of all interacting particles. Besides, 
we keep the main terms of Y^ in ip expansion. The quark zero modes generate the factor similar to 

N f 



Eq. ( |13]) with being changed by the factor ( j^rj an d then in the leading N c order we have 
/ 1 \ Nf f 



dkdl 



Sf 



7o(M) + 



dp 



i(k~l+p)z 



7i(fc,Z) pip) 



ri(k) ^ Ml) ■ 



While providing the Gaussian form for the functional we perform the integration over the auxiliary 
parameter A together with the bosonization resulting in the integration over the auxiliary matrix 
Nf x Nf meson fields [IT 



exp 



, ■ / i J 
Adetl- 



dM exp { i Tr[MJ] - (N f - 1) 



'det[MN c 



As a result the generating functional may be written as 
Z = j^K exp(-iVlnA) J Dp expj- j 



dk nn 



(2k) 



p(-k) [k 2 + M 2 } p(k) 



J DM LtR exp J dz l-(N f - 1) ^ 



'det [M L N c } ^ N f-i + ^det[M R N c ]\ N f 



A 



J D^Dip exp j j ^ f (k)(-k)^ f (k)+i J dz(Tr[M L J + }+Tr[M R r 



(23) 



(2k) 

Now scalar field interacts with the quarks of the different flavours, nevertheless, the dominant con- 
tribution is expected from the tadpole graphs where any pair of the quark fields is taken in the 
condensate approximation as it happens at Nf — 1 

$(*M,(0 - <4(*)Vb(0> = -ir% 9 5(k - I) Tr S(k) . 



1 Q 



As for the condensate itself we obtain it as the nontrivial solution of saddle point equation. For 
example, it is 

(M L , R ) fa = M 5 fg 

for the diagonal meson fields. The dimensionless convenient variables (in addition to Eqs. (0), fllBD ) 
are the following 

M _ 3 / Ap 4 \ jr n J 

~Y P \(2N c p) N f ) 

Then the effective action (Z = J dgdp exp{— V e ff}) in new designations has the form 



V r )n N f- 1 V9N 2 N 2 V r dh 

V eff = N(Nf -l)lng-^(N f - 1)^L_ - j/^ff^) - 2N f N c ^ f ™ \n{k 2 + T 2 (k)} 

and the saddle point equation reads 



np 4 J 7r 4 k 2 + T 2 (A;) np A v np 

The quark current mass in Eq. ([19]) gains the additional factor Nf because the scalar nature of the 
phononlike field requires to match the tadpole quark field condensates of all Nf flavours to every 
vertex 

NfN c p 3 . . 

m fi k > = — =4 C W 7i(«,k) • 

np v 

Table 2 complements the Table 1 with the calculations at Nf = 2 

Table 2. 



D.P. 


M.S.Z. 


p m(o) -t^H) u /; 


p m(o) a /; m y 


4.66- 1(T 3 422 -(427) 3 135 111 


4.61- 1(T 3 416 -(374) 3 118 97.2 12.7 



where f n [MeV] is the pion decay constant and 

f 2 = Ml [ 

U 8 J 



N c Nf rdkT\k)-^T'{k)T{k) + ^{r{k)f f 2 



(k 2 + T 2 (k)) 2 \p c ) ' 

fl [MeV] is its approximated form 

, 2= AyV; rdk T 2 (k) (2_\ 2 
h 8 J tt 4 (k 2 +T 2 (k)) 2 \p c ) ' 

here T'(k) = —^Tr> an d the condensate —i(ip^ijj) is implied for the quarks of every flavour. 

V. Meson excitations of the chiral condensate 

Here we discuss the meson excitations of the chiral condensate adapting the effective Lagrangian to 
the suitable variables. The meson degrees of freedom are introduced into the effective Lagrangian 



fl23|) with Nf = 2 by the following substitutions 



M L = M M L , M L = (l + a + rj) U V , U = exp(z7r a r a ) , V = exp(ia a r a ) 
M R = M M r , M R = (1 + a - n) V £/ f , 



1 A 



where M. denotes the condensate and a, r\ are the scalar and pseudoscalar meson fields, 7r a is the 
isotriplet of the it- mesons, a a is the vector meson isotriplet. The effective action for the meson field 
excitations looks then as 

Veff = - [ % K(p) RAP) * a (-p)+v a (p) RAP) v a {-p) + <r{p) R*(p) <r(-p)+v(p) Rv(p) v(~p)} , 

J 7T 

here R n °., R a a, R a , R v are the inverse propagators of the corresponding particles. Their exact forms 
are listed in Ref. |2j but for the ir- and cr-mesons with the phononlike contributions included they 
are calculated below. 

The contributions generated by the quark determinant if the phononlike fields ignored are well 
investigated in the leading order of the iV c -expansion. In particular, one of the contributions comes 
from the diagram A). It is of the first order in the /i-expansion with one quark loop and one external 




meson fields A4 



Figure 6: Diagram A), 
meson leg in which the meson fields are maintained up to the quadratic terms 

A) 2i Jdzj^ {-Tr S(k)) T(k) [m l h {z) + M*(z) ^} . 

Besides, another contribution comes from the diagram B) being of the second order in /i with one 
quark loop and two external meson legs 




Figure 7: Diagram B). 



B) 2 J *™Tr[S(k)S{t)] T(k, l)T(l, k) [m%(1 - k) l -±^M L gf {k - + (4 7s) - (R, -75)} , 

where we introduced T(k, I) = /r/offc, + 7 — -a 1 h 7i(^> $¥1 7o(&> = G(k)G(l). The determinant 

(np K) ' 

terms of the generating functional at Nf = 2 develop the form 

det M L + det M R = 1 + a 2 + r/ 2 . 

In addition to this diagrams we have the tadpole diagram C) in which the meson fields are maintaned 
up to the quadratic terms and which describes the 1 effect of equilibrium instanton size fluctuations 
when the meson fields are present. These fluctuations influence the IL parameters 0. 

1 K 




Figure 8: Diagram C). 
C) J^Tr (-S(k)) 2 % m f (k) J dz ^M L ff (z) + Mf f (z) ^ 



In particular, we have for the inverse propagator of the 7r-meson field 



r dk r (k, k + p) + T(k)T(k + p) , . ., 

j^M = / ? { - 2 iv e p A ^ _ u^w _ r(t , t + p) r(t + P , k) + 

T 2 (k) 

+2N C Nf L/ + i Nf m f (k) Tr (S(k)) } . 



When the pion momentum goes to zero p — > the first two terms do not generate the massive term 
because of the quark determinant symmetry. As a result, these terms at small values of p originate 
only term proportional to p 2 and we limit ourselves with the approximate result for T = //70 while 
calculating it. Using the expansion obtained in Ref. H] we have 

/d k 
— A m f (k) Tr (S(k)) + 8n p 2 (24) 
7T 

N Nt dk T 2 (k)-kv'{k)T{k) + ^-{T'{k)) 2 /9 ,2 ^ 
with (3n = J p V + r a (fc)) a W ' Eq ' ® combined with the defini " 

tion ( |2"2"| ) can be given as 

From here we have for the 7r-meson mass m 2 = m / ^ji ihUl or jf the relation between 8 and the 

n fin M 

pion constant f n (f 2 = 2 fin) plugged in we obtain 

ml = 2 St y W , (25) 

J 7T 

what displays Gell-Mann-Oakes-Renner relation 

2 (m w + m rf ) i[(u*u) + (c^d)] 



■J 7T 



The factor 2 in the numerator of Eq. fl25|) corresponds just the sum of the u and d quark masses and 
Nf means the summation of condensates. 

The light particle which we introduced and which imitates the scalar glueball properties does not 
affect significantly the SCSB parameters and correctly describes the soft pion excitations of quark 
condensate. Meanwhile, the experimental status of this light scalar glueball is very vague. We believe 
the phononlike excitations could manifest themselves being mixed with the excitations of the quark 



8 The problems related to the effective chiral Lagrangian (when the phononlike fields included) and its symmetries 
will be discussed in the separate publication. 



1 p. 



condensate in the scalar channel. To illuminate the point let us consider the inverse propagator of 
the cr-meson which is given by the contributions of the diagram B) and determinant 



„ , s xr xr f dk (k,k + p) -T(k)T(k + p) , , 4 

«,M = -2n c n s j-, [k ^ r2ik ^ k+ [ J^'[ p)] m k + P ) r(* + p, *) + V . 

Holding the highest terms of the //-expansion only, when T = /ryo, we receive the following result by 
means of the identity (see, f2j) 

R(v) = N N f [dk ink)(k+p), + T(k+p)k,r 
HAP) «c «, J an [k2 + T 2 (mk + p) 2 + T 2 (k + p)] • 

In the course of this exercise we have to include the contribution of the 'shifting' diagram of C) type 
where the field (p (more exactly tp') should be treated as the dynamical one, i.e. 

O 2i» J * Tr (S(k)) ^ / dz {m^z) 1 -^ + M^z) 1 -^ 

Then the term of interacting scalar fields which we are interested in looks like V^ CT = / dz A <p(z) cr(z), 
where A is defined by Eq. At low momenta the diagrams contribute to the effective action (in 

the dimensionless variables) as 

-2 (p 2 + M 2 ) if 2 -2 np a (p 2 + Ml) a 2 + A if a 

with (3 a standing for the kinetic coefficient of the a-meson and M a is its mass. Both quantities 
are derived from the expansion of the inverse propagator in the low momentum region as R a (p) — 
2 njScr (p 2 + M 2 ). Diagonalizing the quadratic form 



-2 {p 2 + M 2 ) V 2 - 2 {p 2 + M 2 a ) a 2 + A <p a 

A 



(presented in more adequate variables a = {nl3 a ) 1 / 2 a, A = ^1/2 ) l ea ds to the following definition 
of the composite particle masses 



M 2 Ml + M 2 ^[Ml-M^ + A? 

M i,2 = g 2 ' ^ 

It is well known that SCSB is inapplicable when we are going to deal with heavy meson masses. 
Specifically, it leads to the wrong predictions for the cr-meson because the corresponding mass ob- 
tained is of an order 1/p ~ 1 GeV (moreover, the adiabaticity approximation is certainly broken 
then). The result Eq. (|26| ) aggravates, in a sense, the situation since the hard component leaves for 
the region of harder masses whereas the light component persists to become lighter. On the other 
hand, the present observations signal rather the existence of the scalar meson of mass about 0.5 GeV 
what apparently does not coincide with the light component M 2 . Nevertheless, our consideration 
having no claims of the quantitative agreement shows the scalar meson could be the mixed parti- 
cle of pretty large width (similar to the superposition of oscillators with the various fundamental 
frequencies) describing the excitations of chiral and gluon condensates. 



VI. Conclusion 

In this paper we have developed the consistent approach to describe the interaction of quarks with 
IL. Theoretically, it is based (and justified) on the particular choice of the configurations saturating 
the functional integral what occurs to be not merely a technical exercise. They are the deformable 



1 V 



(crampled) (anti-)instantons with the variable parameters j(x) and in the concrete treatment of 
this paper we play with the variation of the PP size p(x,z). In a sense, such an ansatz is strongly 
motivated by the form of quark determinant which is solely dependent on the average instanton 
size in the SCSB theory. We have demonstrated that in the long-length wave approximation the 
variational problem of the deformation field optimization turns into the construction of effective 
Lagrangian for the scalar phononlike (p and quark fields with the Yukawa interaction. Physically, 
it allows us to analyze the inverse influence of quarks on the instanton vacuum. We have pointed 
out this influence on the IL parameters as negligible. The modification of the SCSB parameters 
occurs pretty poor as well. In particular, the scale of quark condensate change amounts to a few 
MeV only. Nevertheless, switching on the phononlike excitations of IL leads to several qualitatively 
new and interesting effects. The propagation of the quark condensate disturbances over IL happens 
in this approach to be in close analogy with well-known polaron problem. We imply a necessity to 
take into account the medium feedback while elementary excitations propagating. The intriguing 
conclusion comes with realizing it leads to generation of current quark masses m/ within QCD itself 
with their values entirely corresponding to the conventional phenomenological results for u and d 
quarks. Moreover, the 7r-meson being massless pseudo-Goldstone particle in the standard SCSB 
theory acquires its mass obeying the Gell-Mann-Oakes-Renner relation when the IL deformations 
enter the game. Besides, it hints that fitting the parameters p, n and renormalization constant A all 
together with the alteration of s(p) profile function we might achieve suitable agreement not only in 
the order of magnitude. The difficulties which confronted us here illuminate the fundamental problem 
of gluon field penetration into the vacuum (the instanton vacuum in this particular case) as the most 
principle one. Indeed, it is a real challenge to answer the question about the strong interaction carrier 
in the soft momentum region. Perhaps, the light particle of scalar glueball properties which appears 
inherently in our approach and should manifest itself in the mixture with the excitation of quark 
condensate in scalar channel (cx-meson) is not bad candidate for that role. By the way, it could be 
experimentally observed as a wide resonance. 

Summarizing, we understand our calculation can not pretend to the precise quantitative agree- 
ment with experimental data and see many things to be done. We are planning shortly to consider 
the problem of instanton profile |12| , to make more realistic description of the PP interaction, to push 



our ansatz beyond the long-length wave approximation analyzing more precisely 'instanton Jacobian' 
5 A 
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Appendix 

The contribution of the quark determinant to the IL action is given by the tadpole diagram Eq. (Ell) 
which takes the following form when returned to the dimensional variables (see, Eq. (|15|)) 

nK ) 1/2 .fn\ A /„^.U/2 f ,™(P) f dz P( Z )~Pc 



Aip = A [nK { ip(0) = A (np^K) 1 ' 2 [ dp , . 

p J n J p p 

Then the IL action, Eq. (|3|), acquires the additional term 



(S}= ldzn\(s)-(A> P —^) 



1 Q 



where A' = /i 2 c(/i) and the mean action per one instanton is given by the following functional 
(si) = / dp si(p)n(p)/n with 

si{p) = M + 5 ln(Ap) - \nP 2N < + ^pW - A' (p - p c )/p . 
In order to evaluate the equilibrium parameters of IL we treat the maximum principle 

(e- s ) > (e- So )e-^ s ^ 

adapting it to the simplest version (when the approximating functional is trivial Sq = 0). In a 
sense, this choice of the approximating functional should be a little worse than in Ref. ||. Its 
only advantage comes from the possibility to get the explicit formulae for the IL parameters in lieu 
of solving the complicated transcendental equation. In equilibrium the instanton size distribution 
function n(p) should be dependent on the IL action only, i.e. n(p) = Ce~ s ^ where C is a certain 
constant. This argument corresponds to the maximum principle of Ref. ||. Indeed, if one is going to 
approach the functional @ as a local form (s) = J dp s\{p)n{p)/n where S\{p) = (3{p) + 51n(Ap) — 
\n/3 2Nc + /3£ 2 p 2 np 2 , it makes the approach self-consistent. The functional difference (s) — (si) = 
J dp {s(p) — s\{p)}e~ s ^ /n being varied over s(p) leads then to the result s(p) = si(p) + const 
keeping into the mind an arbitrary normalization. The maximum principle results in getting the 
mean action per one instanton as the IL parameters function, for instance, average instanton size p. 
The corrections generated by the 'shifting' terms turn out to be small and we consider them in the 
linear approximation in the deviation A. The following schematic expansion exhibits how the major 
contribution appears 



(( S + 6) e— s ) _ (s e~') + (6 e~') (s e~ s ) (6 Q - (s 6 e~ s ) (Q 
(Sl) = <e— '> " (e^) + ' (27) 



here 5(A) stands for a certain small 'shifting' contribution and s is the action generated by the gluon 
component only. The last term in Eq. (27) is small comparing to the first one and we neglect it. 
Then it is clear that evaluating the mean action per one instanton is permissible to hold the gluon 
condensate contribution s only (without the 'shifting' term S) in the exponential. Hence we have for 
the mean action per one instanton (si) = J dp si(p)n (p) /no, and no(p) is the distribution function 
which does not include the 'shifting' term 0. It is possible to obtain for the average squared instanton 
size and the IL density that 



n = CC N J 2N ^, (29) 



The 'shifting' term changes the mass of phononlike excitation insignificantly. The equilibrium instanton size 

/ 1 \l/2 

= is equal then to p c — (a + A' p) p, a = ( 1 — J , (3 = 

P=Pc 



as dictated by the condition —g^- 



2jL j\ — a ^iy2p|/^ j and the second derivative of action in the equilibrium point equals to s (p c ) — = + 
= A' 1 V?^ — tt^ — \ ■ Another source of corrections to the kinetic coefficient appears while one considers 



(J oix Z i 

the instanton profile change A —> A + a where the field of correction is a ~ fix — 



This mode could ap- 



pear within the superposition ansatz Eq. (Q) and leads to the modifications of quark zero mode [D{A + a)ij) = 0). 
Fortunately, both corrections to the kinetic term occurs to be numerically small. 

10 In order not to overload the formulae with the factors making the results dimensionless, which are proportional 
to the powers of A, we drop them out hoping it does not lead to the misunderstandings. 
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where the parameter r 2 equals to 
Expanding In p = In p + 



r 2 = P^ 2 np 2 



(30) 



^ — ^ \ p 2 ^ — ^ " ' anc ^ us ^ n S Eq. (p8|) we can show that 



P 



J dp n (p) hip 
J dp n (p) 



lnp + $1(1/) 



/jg n (p)p 
/dp n (p) 



p + $ 2 (i/) 



where $1, $2 are the certain function of z/ independent of p. Besides, the average squared instanton 
size within the precision accepted obeys the equality r 2 p 2 = $(z/), and $(^) is the function of v only. 
Then the mean action per one instanton looks like 

(si) = -2jV c ln/? + (2z/-l) \np + F(u) 

{F(v) is again the function of v only and its explicit form is unessential for us here). Calculating its 
maximum in p we receive 



p = exp 



2N r 



2bN r 



InC 



2u- 1) ' ' 2^-1 
From Eqs. (p8|) , (|30D we find the IL density to be as 

,£^+1/2)1 
2z/ 3 / 2 I» J ' 

and handling Eq. ( p9| ) we determine the constant C. 

The IL parameters occur to be close to the parameter values of the Diakonov-Petrov approach 
[0 and are shown in the following 

Table 3. 




D.P. 


M.S.Z. 


N f pA n/A A (3 


Nf pA n/A 4 (3 


0.37 0.44 17.48 

1 0.30 0.81 18.86 

2 0.24 1.59 20.12 


0.37 0.44 (0.49) 17.48 

1 0.33 0.63 (0.71) 18.11 

2 0.28 1.03 (1.17) 18.91 



Here Nf is the number of flavours, A^ c = 3, and the IL density at A 7^ when the iteration process 
completed is shown in the parenthesis. It is curious to notice the quark influence on the IL equilibrium 
state provokes the increase of the IL density. 

In Table 4 we demonstrate the mass gap magnitude M and the wave length in the 'temporal' 
direction A4 = M~ l . To make it more indicative we show also the average distance between PPs from 
which is clear, in fact, that the adiabatic approximation A>L~^>pis valid for the long- length 
wave excitations of the 7r-meson type. All the parameters are taken at k = 4/5 but the primed ones 
correspond to the kinetic term value of k = 6/3. 

Table 4. 



N f 


MA- 1 


AA 


M'A- 1 


A'A 


RA 





1.21 


0.83 


0.99 


1.01 


1.23 (1.2) 


1 


1.34 


0.75 


1.09 


0.91 


1.12 (1.09) 


2 


1.45 


0.69 


1.18 


0.84 


0.99 (0.96) 



Here the parameters in the parenthesis designate the distance between PP (R 
iteration process is completed. 



n 



-1 / 4 ) when the 



on 
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